Anisotropic Local Stress and Particle Hopping in a Deeply Supercooled Liquid 
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The origin of the microscopic motions that lead to stress relaxation in deeply supercooled liquid 
remains unclear. We show that in such a liquid the stress relaxation is locally anisotropic which can 
serve as the driving force for the hopping of the system on its free energy surface. However, not 
all hopping are equally effective in relaxing the local stress, suggesting that diffusion can decouple 
from viscosity even at local level. On the other hand, orientational relaxation is found to be always 
coupled to stress relaxation. 



Dynamics of supercooled liquid show many fascinating 
properties, namely the non-exponentiahty in various re- 
laxation functions of dynamical variables, like the stress, 
density, composition, polarization and orientation^, |2, 
Hi ^ 1^ ' One often observes a very slow power law de- 
cay in the intermediate to long times Although these 
have drawn attention of experimentalists and theoreti- 
cians and have been extensively studied, many aspects 
still remain ill-understood. It is believed that the nature 
of the microscopic motion undergoes a drastic change 
at a temperature (a temperature substantially be- 
low the melting temperature, Tm)- At the temperature 
Tb, the continuous hydrodynamic type motion, which is 
prevalent at higher temperature changes to discontinu- 
ous motion. It is also believed that is close to the 
temperature where effects of the underlying free energy 
landscape on the dynamics are felt for the first timeQ. 
Tb is also found to be close to the mode coupling temper- 
ature, Tc. This change in the nature of the microscopic 
motion is believed to be the origin of the experimentally 
observed a-j3 bifurcation and also the cross-over be- 
tween the rotational and translational relaxation times 



Due to the complexity of the problem, computer sim- 
ulations have played a key role in augmenting our under- 
standing in this area. In particular, simulations allow one 
to directly look at the microscopic events. The computer 
simulation studies of the stress auto correlation function 
in the supercooled liquid could successfully reproduce the 
power law behavior of the stress auto-correlation function 
(SAF) [1,3. However, in the deeply supercooled liquid, 
one finds that within the simulation time, the relaxation, 
after an initial decay (typically 10-20%) becomes fully 
arrested. The microscopic origin of the subsequent de- 
cay is unclear. The computer simulation studies further 
show that the orientational and translational hopping of 
particles are the only mode present and hence the stress 
relaxation can happen only via hopping. However, since 
the relaxation time is much much longer (could be of 
the order of ms or sec), the computer simulation results, 
(which can explore mostly upto nano second regime and 
sometimes micro second regime), cannot include the ef- 
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fects of these hopping. Therefore, one cannot explore the 
relationship between hopping and the total stress relax- 
ation. 

The experiments on the other hand are successful in 
showing the decay of the SAF. However, except a re- 
cent work using single molecule spectroscopy these 
experiments are macroscopic and do not provide enough 
information of the microscopic motions in the system. 

Majumdar [T^ had earlier discussed the importance of 
local relaxation modes (of wavelength less than the short 
range order) in giving rise to non-exponentiality in the 
stress relaxation function. This work discussed relaxation 
in terms of relaxation within small regions, surfaces and 
also volumes, with progressive lengthening of time scales. 
However, in that analysis the basic mechanism of relax- 
ation was still assumed to be continuous. 

In this Letter we demonstrate for the first time that 
in the deeply supercooled liquid (where hopping is the 
only surviving large amplitude motion), there is a close 
relationship between the local stress and the orientational 
and translational hopping. The local SAF is anisotropic 
and is found to change drastically during the hopping, 
thus showing that the local stress and the hopping of a 
particle are intimately connected. The anisotropy in the 
local stress could be the driving force for hopping. As the 
free energy of the system can be expressed in terms of 
the position dependent stress in a generalized Ginzburg- 
Landau formulation (IL] . the change of the anisotropic 
stress due to hopping should be regarded as the driving 
force for the transitions of the system between different 
minima of the free energy surface. However, not all hop- 
pings are effective in relaxing the stress. 

Our solvent is represented by binary Lennard- Jones 
mixture, which has been extensively studied P, 0, 0, 
IT^ and is known to be a good glass former, and our 
solute probes are prolate ellipsoids. Pressure is kept con- 
stant by Andersen's piston method while in the case 
of temperature, a damped oscillator method has been 
adopted which keeps temperature constant at each and 
every time step |l6|| . The piston mass involved here is 
0.0027(mA/tT^) which is regarded as optimum [TgI It^. 
The interaction between two ellipsoids with arbitrary ori- 
entations is assumed to be given by the Gay-Berne (GB) 
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potential [l^ . 
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where iii 112 are the axial vectors of the ellipsoids 1 and 2. 
f is the vector along the intermolecular vector r = r2 — ri, 
where ri and r2 denote the centers of mass of ellipsoids 
1 and 2 respectively. cr{f,ui,U2) and e(f, ui,U2) are the 
orientation-dependent range and strength parameters re- 
spectively, a and e depend on the aspect ratio k. The 
minor axis of the ellipsoid is equal to the diameter of 
the larger solvent and the major axis is 3 times that of 
the minor axis. Finally, the interaction between a sphere 
and an ellipsoid is accounted for by a modified GB-LJ 
potential given below 
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where 'E' denotes the ellipsoids and 'i' can be 'A' or 'B'. 
The expression for (j{9)Ei is available [l9|. 

The ellipsoid in binary mixture system with the above 
mentioned potential is a well behaved system and it 
can also exhibit the experimentally observed anoma- 
lous viscosity dependence of the orientational correla- 
tion time [13 ■ Four ellipsoids were placed far from each 
other in a binary mixture of 500 particles with num- 
ber of 'A' particles, Na = 400 and number of 'B' type 
particles Nb — 100. The reduced temperature is ex- 
pressed as, T*{—kBT/eA), the reduced pressure as, P*{— 
Po'a/^aa)- and the reduced density as p*{=pa\). The 
time is scaled by r = {mA<j\/ eAA)- The time step of 
the simulation is .002 r and the system is equilibrated 
for 1.5 X 10^ steps and the data collection step is 5 x 
10^. The studies have been performed at T*=0.8 and the 
P*=6 and 10. 

At P* =6, both hopping and continuous motions ex- 
ist in the system, thus the stress relaxation could not be 
directly correlated with the hopping. At P*=10, only 
microscopic motion that survives is hopping. In a recent 
study, we have reported observation of correlated trans- 
lational and orientational hopping at this pressure. 
After extensive simulations we could find only two dif- 
ferent kinds of motions. The translational hopping was 
either associated with correlated hopping of 5-6 nearest 
neighbors or it exhibited a motion in a ring like tunnel. 
While it is possible that other types of motion like iso- 
lated hopping can exist, we could not find them. The 
hopping rate was found to be 2 x 10^ where both the 
type of motions occurred with almost equal frequency. 
In the following we focus on the stress relaxation and its 
relation vis-a-vis hopping at P* = 10 and T*=0.8. The 
reduced density of the system is 1.41. 

In figure 1 we show two different kinds of spatial hop- 
ping observed in our simulations. Both the hoppings 
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FIG. 1: (a) Shows the displacement of the 1st tagged elhpsoid 
over 1000 r. There is a hopping of the ellipsoid around 300 
T. In the inset we plot the orientational correlation function 
obtained at different intervals, (b) Same as figure 1 (a) but 
for the 2nd tagged elhpsoid. The plots are at P* = 10 and 
T*=0.8. 



are associated with orientational hopping. In the inset 
we also plot the orientational time correlation functions 
(OCF), before, during and after the hopping. Figure 
1(a) shows the trajectory of the first tagged ellipsoid and 
inset shows its orientational time correlation function, 
< P2{u^{Q)u^{t)) > /< P2{u^{0)u^{0)) >, were P2 is the 
second order Legendre polynomial. The hopping takes 
place in 20t and the displacement is 0.5cr. Here the el- 
lipsoid hopping is accompanied by hopping of 5-6 of its 
nearest neighbors. The OCF decays only during the pe- 
riod of hopping. In figure 1 (b) the trajectory of the 
second tagged ellipsoid and in the inset its orientational 
time correlation functions are shown. Note that in this 
case the displacement of the particle is large (l.lcr) and 
it also takes place over a very long period of time (50 r). 
Here we find that the tagged particle is moving in a ring 
like tunnel. Although the orientational correlation func- 
tion decays during the hopping, it's decay is less when 
compared to that of the 1st tagged particle. 

The motions in a string like and ring like tunnel in a 
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FIG. 2: (a) The different components of the local stress auto- 
correlation (obtained from the sum of the stress of the 1st 
ellipsoid and its nearest neighbors) function obtained before 
hopping, between 320-330 r. (b)is the same as 2 (a) but the 
components of the SAF are obtained after hopping between 
(340-350r). The plots are at P*=10 and T*=0.8. 



FIG. 3: (a)The different components of the local stress auto- 
correlation (obtained from the sum of the stress of the 2nd 
ellipsoid and its nearest neighbors) function obtained before 
hopping, between 1590-1630 t. (b) is the same as 3 (a) but 
the components of the SAF are obtained after the hopping, 
between 1670-1710t. The plots are at P*=10 and T*=0.8. 



deeply supercooled liquid have been reported earlier by 
other authors, although they did not discuss this large 
disp lacements over a long time for the ring like motion 

In the following part we discuss in details the local 
stress relaxation associated with these two different kinds 
of hoppings. The local stress around the ellipsoid is 
obtained by summing over the stress on the ellipsoid 
and its nearest neighbors. The stress has six compo- 
nents and the stress auto-correlation function is given 
by < aij{0)aij{t) >, where i,j=x,y,z. In case of only 
Lennard- Jones fluid, o-y = ajt, but for particles interact- 
ing via GB and GB-LJ potential this is not so. 

Figures 2 (a) and (b) show the SAF around the 1st 
tagged ellipsoid, before and after the period it is hopping, 
respectively. Before the hopping there is an anisotropy 
of the stress. The xy and yx components of the stress are 
much larger than the others and also their correlations 
cease to decay. This anisotropy leads to a hopping of the 
ellipsoid, mostly in the z direction. During the hopping 
there is a relaxation of the SAF in xy and yx components 
and after hopping all the components relax. Note that 
the total stress in all the components are also lower. 

Figures 3 (a) and (b) show the SAF around the 2nd 
tagged ellipsoid, before and after the period it is hopping. 



respectively. Before the hopping there is an anisotropy 
of the stress. The yz and zy components of the stress 
are much larger than others, and also their correlations 
cease to decay. This leads to the hopping of the ellipsoid, 
mostly in the x direction. During and after the hopping 
there is an exchange of the stress. After the hopping 
although the yz component of the stress relaxes and also 
the t=0 value of all the components reduces, the SAF 
in the xz and zx components cease to decay. Thus, this 
kind of motion in a ring like tunnel does not lead to the 
relaxation of all the components of the stress. 

The orientational correlation function of the 1st tagged 
particle relaxes more ( inset of figure 1 (a)) compared to 
that of the 2nd tagged particle (inset of figure 1(b)) when 
computed in the respective intervals where they are hop- 
ping. From figure 2 and 3 we found that the local stress 
relaxation takes place when the 1st ellipsoid is hopping 
where as when the second ellipsoid is hopping although 
there is an exchange of stress between it's components, 
the local SAF does not completely relax. There is a direct 
connection between the local stress and the orientational 
relaxation functions, implying that rotation and viscosity 
are coupled even in a localized region. 

In order to understand what happens to the surround- 
ing of the local region when there is a relaxation of stress 
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due to hopping wc have studied the stress auto correla- 
tion function of a bigger region of 2 cr around the 1st 
elhpsoid.We found there are about 62-67 particles in this 
region. Although there is an anisotropy of the compo- 
nents of the stress in this larger region, this anisotropy 
cannot be correlated with the direction of hopping. Some 
of the components of SAF build up immediately after 
hopping and in a later time it relaxes. Similar analysis 
when done in a even bigger region shows that it takes 
longer for the stress of that region to relax and also the 
effect of the hopping is less. 

In conclusion we demonstrated that the direction of 
the hopping of the tagged particle is determined by the 
anisotropy in the stress. Anisotropic stress relaxation is 
different when there is a many-particle hopping and there 
is a motion in a ring like tunnel. Although there is an 
exchange of stress between the components due to the 
particle motion, the stress relax;ation is less in a ring like 
motion. Interestingly, the effect of hopping is found to 
spread over the adjoining region like ripples with the am- 
plitude decreasing with increasing distance from the el- 
lipsoid. We note that in the case of the second tagged el- 



lipsoid (Fig. lb) although it translates more, the stress re- 
laxation during its hopping is less. Thus suggesting that 
translational motion and viscosity are decoupled even in 
a localized region. On the other hand, the orientational 
relaxation and also the stress relaxation is more for the 
first ellipsoid. Thus suggesting that the orientational mo- 
tion always remains coupled to viscosity. This is in agree- 
ment with the experimental results and in fact provides a 
microscopic explanation of the results which are known 
for a long time. There can be an apparent connection 
between the stress tensor and the momentum circula- 
tion. Thus it is possible that the non-decaying SAF im- 
plies that momentum circulation exists in a deeply super- 
cooled liquid Since the anisotropic stress contributes to 
the free energy of the system, a change in the anisotropy 
drives the system from one free energy minimum to the 
other. When the anisotropy in the stress disappears and 
all the SAF relaxes then the system definitely moves to 
a lower free energy minimum. 
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